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We propose a general method to find exact ground state solutions to the SU(1V) invariant matrix 
model arising from the quantization of the 11-dimensional supermembrane action in the light-cone 
gauge. We illustrate the method by applying it to lower dimensional models and for the SU(2) group. 
This approach can be used to find ground state solutions to the complete 9-dimensional model and 
for any SU(1V) group. The supercharges and the constraints related to the SU(2) symmetry are 
the relevant operators and they generate a multicomponent wave function. In the procedure, the 
fermionic degrees of freedom are represented by means of Dirac-like gamma matrices. We exhibit a 
relation between these finite N matrix theory ground state solutions and SUSY quantum cosmology 
wave functions giving a possible physical significance to the theory even for finite N 
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I. INTRODUCTION 

As is known, M-theory should be an eleven dimen¬ 
sional supersymmetric quantum theory from which all 
superstring theories can be deduced and it can be con¬ 
nected also to 11-dimensional supergravity. This idea was 
born since it was realized that the five consistent theo¬ 
ries of strings and 11 -d supergravity are related among 
them through several types of dualities ElBj ■ There are 
interesting approaches and calculations that provide the 
possibility to understand some important features of this 
fundamental theory. There is a concrete formulation of 
M-theory which has been claimed has the appropriate 
expected properties and it is based on the matrix model 
conjectured in (BFSS 0)- This matrix model is a super- 
symmetric quantum mechanics where the bosonic degrees 
of freedom are a finite set of N x TV Hermitian matrices 
Q. We can associate a time dependent matrix to every 
point in space and these matricial coordinates are obvi¬ 
ously noncommuting. The existence of dynamical objects 
of dimensionalities higher than one, like Dp-branes, has 
left aside the fundamental status of particles and strings 
and the origin of this matrix model and its connection to 
M-theory arises from the dynamics of DO branes mi¬ 
lt was conjectured in [l| that M-theory is the TV — y oo 
limit of the maximally supersymmetric matrix Hamilto¬ 
nian emerging from the DO-brane matrix model. The 
conjecture was extended in 0 showing that matrix the¬ 
ory is also meaningful for finite N, so the solutions to the 
matrix model for N finite are relevant. 

In the recent years, research has been done on another 
matrix model arising from noncommutative gauge the¬ 
ories El , it is called the IKKT matrix model and it 
can also be related to a large N reduced model of 10- 
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dimensional non-commutative Super Yang-Mills theory. 
In the IKKT framework, the parameter 6 that measures 
the non-commutativity between the coordinates of space- 
time becomes dynamical and in the semiclassical limit it 
defines a Poisson structure such that the metric of space- 
time itself emerges naturally jl^4l3 | . This matrix model 
can be seen as a N = 4 noncommutative Super Yang- 
Mills theory on R|. The path integral quantization of 
this matrix model is naturally defined by integrating over 
the space of matrices and in a semiclassical way, it repro¬ 
duces some interesting features of curved spacetime |l5[ . 
So, it has been claimed is a good candidate to be a quan¬ 
tum theory of g ravity together with other fundamental 
interactions [161 ] . The presence of branes in this theory 
El becomes also important as it is in the matrix theory 
of BFSS. The quantization approaches are different for 
the two matrix models briefly described above, but they 
have the same common Dp-brane world origin and both 
involve the most important characteristic features of a 
unified theory of physical interactions including gravity. 

Another, and somehow independent, matrix model 
was encountered in the quantization of the classical 11- 
dimensional supermembrane action in the light cone 
gauge (LCG) [18l - l2l| . These matrix models were stud¬ 
ied in the context of supersymmetric quantum mechanics 
years before the quantization of the supermembrane was 
finally achieved [22|. The action of the supermembrane 
is the generalization of the Green-Schwarz action for the 
superstring for a membrane moving in a d-dimensional 
target superspace, and when d = 11 this can be directly 
related to 11-dimensional supergravity [23|. It is possible 
to regularize the theory because of the existence of a finite 
dimensional Hilbert space where the area preserving dif- 
feomorphisms (APD) invariant states can be expressed. 
This APD invariance is translated into a SU(iV) sym¬ 
metry in the regularized theory which is equivalent to a 
Super Yang-Mills quantum mechanics matrix model 0 
l2p| . The Hamiltonian and the associated ground state 
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of these matrix models of supersymmetric quantum me¬ 
chanics with SU(7V) invariance were studied in [22] and 
the matrix model arising from the supermembrane action 
is exactly one of those. In this type of supersymmetric 
models it is of great importance to find the ground state 
solution [2l| 122, l2J, 125| , and for the matrix model we have 
been briefly mentioned here the existence of such state, 
even for N finite, has been found only perturbatively. 

These two matrix models, the one related to DO-branes 
and the other related to the supermembrane, are not in¬ 
dependent. There is a connection between membranes 
and DO-branes first speculated in [26] where it was con¬ 
sidered that membranes could be regarded as collective 
excitations of DO-branes. This connection was treated 
in great detail in jlj. The noncommutative space defined 
by the phase space coordinates where the membranes ex¬ 
ist has a basic minimum area, this quantum unit cells of 
space are precisely the DO-branes, each one with a lon¬ 
gitudinal momentum of 1/iV. The longitudinal direction 
is the one compactified in the LCG description. This 
relation of duality connects the two matrix models. 

The first purpose of this work is to exemplify the gen¬ 
eral method by finding solutions to the ground state of 
matrix theory reduced models arising from the quantiza¬ 
tion of the classical supermembrane action. Our proce¬ 
dure is based on; i) We consider finite N matrix mod¬ 
els, ii) The relevant operators are the supercharges and 
those related with the SU(iV) symmetry. It is no more 
necessary to consider the Hamiltonian operator because 
through the canonical algebra it is a consequence of these 
constraints, iii) The fermionic degrees of freedom are rep¬ 
resented utilizing a Dirac-like gamma matrix representa¬ 
tion ( 27 ], Ell • By means of these, the associated super¬ 
charges and the operators related to the SU(iV) symme¬ 
try become matrix differential operators acting on a mul¬ 
ticomponent wave function. One generates a Dirac-like 
quantum matrix model formalism. 

With our proposal one can solve any finite N matrix 
model and for all the bosonic and fermionic degrees of 
freedom. Even though finite N matrix models are phys¬ 
ically relevant 0, it is not clear which particular N 
should be considered. For large N, our method will gen¬ 
erate a large number of first order coupled differential 
equations. In order to illustrate how to apply our pro¬ 
posal, we restrict ourselves to reduced models, for the 
SU(IV = 2) invariance group. Besides, for the sake of 
simplicity and as a way to demonstrate how this pro¬ 
cedure works, we will make further assumptions in the 
space of bosonic variables. We will then freeze out some 
bosonic degrees of freedom. By doing this one cannot 
assure that this reduced models will capture all features 
of the quantization of finite N matrix theory (“for a re¬ 
lated discussion in the context of quantum cosmology 
see; 0”). Some of the models studied can however be 
exactly solved illustrating us about the kind of ground 
state solutions to be expected for a larger N model. The 
method we propose can be directly extended to higher di¬ 
mensional models and for any SU(iV) group. This kind of 


representation for the fermionic variables has been used 
in a first and several models of SUSY quantum cosmology 
0, EM 113113 leading to pretty much interesting features 
regarding every physical supersymmetric system treated 
in this way. The fermionic degrees of freedom can also 
be represented by differential operators [H|-|37[. A rel¬ 
evant point is that it is possible to solve the complete 
operators system containing both, the bosonic and the 
fermionic sectors. The fermionic information will then 
be encoded in each one of the entries of the multicompo¬ 
nent wave function and the physical information regard¬ 
ing the variables of the system will be given by means of 
the components of this wave function. 

It has been argued that M-theory is a promising can¬ 
didate to be a quantum theory of gravity, it is then ex¬ 
pected that the matrix models related with it, contain 
quantum and classical gravitational information. The 
same authors of the conjecture have shown that a spe¬ 
cific compactification of matrix theory correctly describes 
the properties of a Schwarzschild black hole |38l - l41j and 
other works related to matrix theory black holes and its 
thermodynamic properties are |42 h45| . 

It has also been claimed that classical cosmological 
models can be deduced directly from the matrix theory of 
BFSS 0, H3] and specifically in [13 , using a homothetic 
ansatz, the authors were able to derive the Friedmann 
equations from the bosonic sector of the classical Hamil¬ 
tonian of matrix theory. Following this kind of ideas, the 
interesting point to us here is that we show that some su¬ 
persymmetric cosmological models are then hidden in the 
SU(./V) invariant supersymmetric matrix model as well. 
It is then expected to find supersymmetric quantum solu¬ 
tions that can be interpreted and related to SUSY quan¬ 
tum cosmology models. This is a particularly interest¬ 
ing byproduct of our reduced matrix models and shows 
that these reduced models are also supersymmetric given 
that their corresponding SUSY cosmological models are 
supersymmetric by construction [3ll f32j j . Also recently 
[48j , to describe the quantum dynamics of the supersym¬ 
metric Bianchi IX cosmological model 0] it has been 
proposed to fix from the start the six degrees of freedom 
describing local Lorentz rotations of the tetrad. The op¬ 
erational content of the supercharges and Hamiltonian 
revealed a hidden hyperbolic Kac-Moody structure which 
seems to support the conjecture about a correspondence 
between supergravity and the dynamics of a spinning par¬ 
ticle moving on an infinite dimensional coset space. At 
the same time we show that some reduced matrix models 
correspond to exact SUSY quantum cosmology models 
that arise from supergravity. All these relations point 
out to a web of connections between these different top¬ 
ics and make it of interest to search for a method to find 
ground state solutions in finite N matrix theory mod¬ 
els. Here, we present a general procedure to construct 
them. We will also be interested in relations between 
the solutions of these ground state matrix theory models 
and wave functions corresponding to SUSY quantum cos¬ 
mology. One of the wave function solutions obtained in 
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our reduced models resembles a wave function solution 
obtained in the superfield approach to supersymmetric 
quantum cosmology and it is then argued that the ma¬ 
trix model we consider provides the relevant physical in¬ 
formation, as classical cosmology seems also to arise from 
bosonic matrix models !46, 47]. 

The work is organized as follows. First in section II we 
show the matrix model and its description as a supersym¬ 
metric SU (N) invariant quantum mechanical model. In 
section III we illustrate the way in which we look for solu¬ 
tions with the matricial representation for the fermionic 
degrees of freedom in a 2 -dimensional reduced model and 
for the SU(2) group. In section IV we show the solutions 
we obtain for some particular assumptions on the bosonic 
variables. As we pointed out, this method could be di¬ 
rectly extended to the complete 9-dimensional theory and 
for any SU(N) group, although for large N powerful com¬ 
putational tools would be needed. Then in section V we 
make an extension to a 4 dimensional model. It is inter¬ 
esting to remark that some of the 2 -dimensional solutions 
are also solutions to the 4-dimensional model. In section 
VI we show first that in this extension we encounter, 
exactly, a model guessed in [ 2 d] and used to analyze 
the spectrum of the matrix Hamiltonian of the quantum 
supermembrane action and then we exhibit some solu¬ 
tions to this extended model. In section VII we discuss 
a connection of one of the solutions to the extended 4- 
dimensional model with SUSY quantum cosmology. Sec¬ 
tion VIII is devoted to conclusions. 


II. MATRIX MODEL 

As we mentioned in the introduction, the model arising 
from the quantization of the classical supermembrane ac¬ 
tion, once the regularization is achieved, becomes a model 
of supersymmetric quantum mechanics [lj that can also 
be obtained from a dimensional reduction of a super 
Yang-Mills theory with gauge group SU(iV) [TMIHl, the 
Hamiltonian of this matrix model [2l| is given by 


1 

2' 

igH r 

2 Jabc^-a 


.D^CAcft 


to, n = { 1,2,..., 9} , a, &, c = { 1,2,..., N 2 
a,/3 = { 1,2...,A f} , 


(1) 


1 }, 


where and i r™ are the bosonic degrees of freedom 
and its related momenta respectively, and are the 
components of the SU(IV) valued matrices 4> m given by 
<F m = </>™r a where r a are the basis group elements. Then 
a is a group index and m is a dimension index. A aa 
represent the fermonic degrees of freedom and T m are m 
Dirac matrices obeying the Clifford algebra { T m , T n } = 
25 mn . The SU(iV) group structure constants are f abc . 
The algebra between the supercharges as function of the 


total Hamiltonian and the supercharges themselves are 
given by 


{Q a ,Q f 3 } = 25 a pH + 2g{T n ) a ^ n a G a , ( 2 ) 

Q a = (T m A a ) Q <* + igf abc { E"A a ) 0 CC, 

where E mn = — |[r m ,r n ]. In the anticommutator of the 
supercharges it appears the operator (constraint) related 
to the SU(AT) invariance G a which in terms of the bosonic 
and fermionic variables is given by 


Ga = fabci^T - yA bQ A ca ). (3) 

Now, we want to realize the quantization. And as we 
are working with a supersymmetric quantum mechanics 
where the Hamiltonian is the square of the supercharges 
and it is positive semi-definite, the state |T) which obeys 
-ffl’F) = 0 is automatically the ground state wave func¬ 
tion. Also, the ground state should be a gauge invariant 
state, it should satisfy G a | v I') = 0. The commutator and 
anticommutator relations for the bosonic and fermionic 
degrees of freedom will be 


[(t>™,Ttb] = ih5 ab 8mn, { A aa , A b0 } = 6 ab 6 at 3 . (4) 

The supercharges arise as the square root of the Hamil¬ 
tonian, and as it can be seen in the algebra eq. 0 , we can 
solve the equations <5 Q |4/) = 0 together with G a |'F) = 0. 
Then the ground state satisfying this, will respect all the 
symmetries; gauge and exact supersymmetry. The rep¬ 
resentation of the canonical momenta corresponding to 
is, as usual, 7 r™ = — ih and given the algebra 
for the fermionic variables eq. 0 . we have the freedom 
to represent the A aa as matrices. This kind of represen¬ 
tation was used first in SUSY quantum cosmology ETa . 
This is an essential tool [27[ l28l i30i. 0 in this work and 
as we will show, we will be able to find exact solutions 
to some reduced lower dimensional models. The method 
we will present can be applied to search for solutions to 
the 9-dimensional supersymmetric matrix model and for 
any SU(iV) group. However, as it will be seen, for these 
dimensions and a large N this is a difficult task because 
one would need to solve an extremely large number of 
coupled partial differential equations. 


III. THE METHOD OF SOLUTION 

To begin with, we will try to find solutions to the 
ground state by solving the corresponding equations 
Qa | tk) = 0 and G a | 4/) = 0 (the Hamiltonian constraint 
is then satisfied) for a SU(2) symmetry; a = 1,2,3 and 
for a 2-dimensional reduced model, to = 1,2. Hence, 
we will have six <fi™ bosonic degrees of freedom and the 
same number for the fermionic ones. In the variables 
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A aa , a is a group index and the a index represents the 
number of supercharges which is in accordance with the 
number of components of a spinor in d-dimensions. Then 
with d = 2 , we have two supercharges and consequently 
six fermionic variables. The dimension M of the square 
M x M matrix representation A aa depends on the num¬ 
ber of fermionic degrees of freedom. So, in this first toy 
model M = 8 . The representation we use for the A aa 
matrices that satisfies eq. © is the following 


Q1 — Ah7Tj + A2l7r| + A3i7r| (8) 

+ Ai 2 7 Ti + A 2 27 T 2 + A327I-3 

-dAl2(</>2^3 - 4>\<t>l) 

- gA 22 {(j)\(t>\ - <t>\4>l) 

- gA 32 0 ^! - , 


^22 


0 ~ 

<1 

11 

A -^ 

A21 = 

±! ,n 

71 ’ 

. 3 

= 71 ’ 

Asl - ^ 

A32 = 

a a 3 

-p a ' 5 . 

72 

= 7^ ® a 1 , 

A 4 = I (g) ia 2 

, A 5 = 

/ ® ia 3 


(5) 


/u = 0,1,2,3, 


Q 2 = AllTrj + A 2 l7T2 + A 3 l7T3 

— Ai27T^ — A 2 2^2 ~ ^ 32^3 

+ 5All(^>2^3 - 4>\<t>2) 

+ 5 A 2 l(</> 30 i - <j)\l tf) 

+ 5^-31 (</>! 02 - , 


where a 1 are the Pauli Matrices, 7 ^ are the Dirac ma¬ 
trices in a Majorana representation and I is the 4x4 
identity matrix. The algebra in eq. © is not affected 
if we choose the plus, or minus, sign on every matrix in 
eq. ©• Notice that we have three kinds of indices, those 
corresponding to the supercharges a = 1 , 2 , to the group 
a = 1,2,3, and the ones associated to the dimensions 
m = 1,2. Even when we have three kinds of indices we 
will denote all of them, as in eq. ©, with the same kind 
of number. As an example, the first of the A aa matrices 
is explicitly 


and the three operators G a (a = 1 , 2 , 3) related with the 
SU(2) gauge symmetry are given by 

Gl =^2 7r 3 + 4>2^3 ~ <\>\A ~ <\>\A (9) 

— ih{ A 21 A 31 + A 22 A 32 ), 

G 2 =4>1'k\ + (fenl - <p\irl - 

— ifi(A 3 iAn + A 32 A 12 ), 
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0 

0 

0 

0 

0 

0 

-* \ 


0 

0 

0 

0 

0 

0 

—i 
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0 

0 

0 

0 
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0 

0 

. ±1 

0 

0 

0 

0 

i 

0 

0 

0 

An “7i 

0 

0 

0 

—i 

0 

0 

0 
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0 

0 

—i 

0 

0 

0 

0 
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0 

i 

0 

0 

0 

0 

0 
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\ * 

0 

0 

0 

0 

0 

0 

0 / 


and for this model we choose the following two T m ma¬ 
trices 


T 1 



1 

0 



(7) 


The two supercharges Q a of eq. © are given by the 
following operators 


G 3 =(j)\nl + (ftirl - (j) 271-3 - 4>1-kI 
— ih( A 11 A 21 + Ai 2 A 22 ). 


Using the matrix representation of eq. all of them 
are linear matrix differential operators acting on an eight 
component wave function IT). We can restrict this model 
even more, setting some of the bosonic variables equal 
to zero, or making some identifications between them. 
These will give specific configurations. Using the anti¬ 
commutator of the supercharges in eq. © and eq. © we 
can see that it is not always necessary to impose all of 
the three G a operators, and that we can solve the prob¬ 
lem for one single operator Q a . For every G a |\I/) = 0, 
and Q a |T) = 0 , we have eight coupled partial differen¬ 
tial equations for the eight components of the wave func¬ 
tion |T) = (^ 1 ,^ 2 ,^ 3 )^ 4 ,^5 ,4>6,4>7,4>8)- For example 
the equation < 5 i|T) = 0 leads to the eight equations 
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[3/i(</>) - ^ 3)^1 + [*tt| - n{ + 3/ 2 (0 )]^ 2 + [^2 - 3/ 3 ((/>)]V’4 - [iA + TrfJV’s = 0 , ( 10 ) 

[ 3 / 2 ( 1 ?) - *tt 3 ~ 7r i]V’i + [tt 3 - 3/i(0)]-02 + [^2 - 3/ 3 ((?)]V’ 3 - [*ttJ + TrfJV’T = 0, 

[A - 3/ 3 (<?)]^2 + [ 3/1 (?) - 71-^3 + [ml + A - 3/ 2 (?)]V>4 + [*7T? + 7i|]V>6 = °> 

[tt 2 - 3/ 3 (?)]^i + [A - i-A - 3/ 2 (?)]V'3 + [A - 3/1 o?)]^ + [*7T? + 7r|]V»5 = 0, 

[772 - mf]^ + [ 3/1 (?) - TrgJV’s + [iA - TTi + 3/ 2 (?)]^> 6 + [A - 3/3(?)]V , 8 = 0 , 

[A - «r?]^3 - [*7rf + A - 3/ 2 (?)]^5 + [tt 3 - 3/i(?)]'06 + [A - 3/ 3 (?)]^r = 0, 

\m\ - A]^2 + [A - 3/ 3 (?)]V’ 6 + [3/i(?) - + [*7r 3 + A ~ = 0, 


[iA - + [A ~ 3/ 3 (?)]V>5 - [iA 


where the functions / a (</), a = 1,2,3, depend on the 
components of the matrices 4> m and give rise to the po¬ 
tentials appearing in the Hamiltonian. These functions 
are given by 

/i(?) = (AA ~ AA), ( 11 ) 

./*2(?) = {AA - AA), 

/3(?) = (A A -AA)- 

We can then make the mentioned simplifications to find 
a model with a particular potential, or even easier, for a 
null potential. This last is the first one for which we will 
find a solution. 

IV. EXACT SOLUTIONS 

Here we will show the explicit form of the ground state 
wave function for a few models using some specific con¬ 
figurations of the components of the <j>™ variables. 

A. Solution 1 

The model is the following. We choose the bosonic 
variables to be (in units where H = 1) 

A =x, A = A = 0, A = y> A = A = o, ( 12 ) 
A = A = 0 j = 2 ’ 3 ’ 

A = A = 0 / = 2 ’ 3 , 

and as a consequence of this choice we have 

/i(?) = 0, (13) 

/ 2 (?) = 0 , 

/ 3 (?) = 0. 


A + 3/ 2 (?)]^7 + [713 - 3/i(?)]V>8 = 0 , 


Hence, in this configuration there will be no bosonic po¬ 
tential and the G a operators take the simpler form 

G\ = —t(A 2 iA 3 i + A 2 2A 32 ), (14) 

Gi = —i(A 3 iAn + A 32 Ai 2 ), 

G 3 = —i(AnA 2 i + A 12 A 22 / 

One should apply all the three G a operators to the wave 
function, however given eq. © and eq. © we have 
to consider only the operator Gi acting on the wave 
function 141), and consistently it is sufficient to apply 
also one of the supercharges, namely Qi|41) = 0. Re¬ 
lations between the eight components of 141) arise and 
the wave function as a row vector has the form |4t) = 
(V’lj — # 2 , — #i, — 4 > 2 , —i>i), where we mean 

by this that, 7/3 = = *^> 2 , and so on. Because 

of the form of the differential equations for this two di¬ 
mensional model, it is easier to solve them using polar 
coordinates x = r cos <p,y = r sin <f>. The exact solution 
for the two independent components of the wave function 
results in 

t/>i = ±ar K e iKC * ± br~ K e~ iK *, (15) 

t/> 2 = Tar K e iK<t> =F br-^e^ . 

Here k is a separation constant and a, b are arbitrary 
constants and we can choose them to have a normalizable 
solution. 


B. Solution 2 

The following model is chosen in such a way that the 
bosonic part of the potential in its corresponding Hamil¬ 
tonian takes the form of a one dimensional harmonic os¬ 
cillator, this is 
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(j)\= x, (j)\ = c = constant, (16) 

02 = 03 = 0i = 03 = o, 

/i( 0 ) = cx, 

/ 2 ( 0 ) = 0 , 

/ 3 ( 0 ) = 0 . 

Following the same procedure as in the previous solution, 
we will have two independent components, 0 i and 02 , for 
the wave function |'F) = (0i, 02 , 02 , 0i, 0i, 02 , 02 ,0i). 
The explicit form of these two components is 


where Ci,C 2 are constants. We see that the constant 
C that relates 0i and 05 depends on the particular 
choice of the constants Ci,C 2 in the linear combination 
of exponentials. The most general solution of 01 can 
have eight terms, each one with a proportional constant 
c,j. j = 1 , 2,..., 8 and the constant C would depend on all 
these Cj. 

D. Solution 4 

This model has three bosonic components, one is con¬ 
stant, with a potential different from zero. These are 


0, = ( a + b)e^ + {a-b)e~^, (17) 

0 2 = i(a + b)e^ — i(a — b)e ~^, 

where a, b are arbitrary constants, if we choose them to 
be a = — 6 , then we get a normalizable wave function. 

C. Solution 3 


01 = X, 02 = 0, 03 = z, (20) 

01 = 0 , 02 = c, 03 = 0 , 

/i( 0 ) = cx , c is a constant, 

02 (0) = -CZ, 

03 (0) = 0. 


The next model has also a null potential, but it has 
three independent bosonic variables 


01 = X; 02 = Zb 03 = (18) 

0 ? = 0 , 0-1 = 0 , 0 | = 0 , 


0i (0) = 0, 

02 ( 0 ) = 0 , 

03 ( 0 ) = 0 . 

We can see from eq. 0 and eq. 0 that for this 
particular configuration it is not necessary to impose 
any G a \ijj) = 0 because 0 2 = 0 , i = 1,2,3. We will 
try to find a solution to Q 2 I 0 ) = 0. This is possi¬ 
ble and the eight independent components of the wave 
function reduce to only two, these are 01 and 0 s, and 
the wave function would take the following form |\I/) = 
(01,01,01,01,05, -05, -05,05). Besides, there is a cou¬ 
pling between these two components, so they are also 
related. The form of the solution of 0i is a linear combi¬ 
nation of the exponential functions that come out from 
every particular choice of signs (±) in the exponential 
factor exp[±\/a "0 ± yfby ± iy/a + b z], where a, b are 
constants and because of the coupling, 05 results pro¬ 
portional to 0i, so 05 = (70i. The form of the constant 
C depends on the particular form of 0i. For instance 0i 
and 05 could be the following 


0 1 _ Cie Vax+Vby+i^/a+bz _|_ C2 g y/ax+y/by-iy/a+bz ^g-j 


05 = 


(ci + c 2 )(Vb- iy/a) 
(ci - c 2 )y/a + b 


01 


After imposing the only non-trivial condition G 2 1 tF) = 0 
and solving the equations Qil’F) = 0 we have a wave 
function solution with just one independent component, 
0i. The form of 0i and its relation with the rest of the 
components is 


01 

01 


= exp 
= 04 = 


i tan 

: 05 = 



08 - 02 


01 = *02 
= 03 = 06 


( 21 ) 


07, 


and a constriction for this model arises, x 2 + z 2 = 

We have seen so far that these reduced lower dimen¬ 
sional models have exact solutions. We also found that 
for some of the configurations we have used, the arbitrary 
constants in the components of the wave function solu¬ 
tions can be chosen so that they are normalizable. We 
could make an extension of the method to a higher di¬ 
mensional model but, as the dimension of the matrix rep¬ 
resentation for A aa grows exponentially with the space 
dimensions, it would be hard to handle the complete 9- 
dimensional problem. It will be needed a powerful com¬ 
putational method to solve the complete matrix model, 
even for the SU(2) group and even harder for a larger 
group. In the next section, we solve a particular 4- 
dimensional model. Our intention is to show that the 
method can be directly extended to more dimensions and 
to any SU(iV) group even though, as we mention, these 
extended models will be difficult to handle. 


V. SU(2) 4-DIMENSIONAL EXTENSION 

In the 4-dimensional extension for the SU(2) group 
there are 12 bosonic variables 0™, to = 1,2, 3,4. a = 
1,2,3, and with d = 4 the number of supercharges is 
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4. The number of fermionic variables is also 12 and the 
dimension M for the matrix representation of these vari¬ 
ables is M = 64. The (64 x 64) representation we choose 
for the twelve matrices A aa that satisfies the algebra eq. 
ED is the following 

An = (il/v^er 1 gilgilglglgl, (22) 

A 12 = (±l/\/2)tr 2 g)lg)lglg)lgl, 

Ai 3 = (±l/y/2)a 3 <S)cr 1 <S)liS)liS)l<S>l, 

Ai4 = {±1/V2)a 1 g)CT 2 glg)lglg)l, 

A 2 1 = (±l/\/2)cr 3 <g) cr 3 g CT 1 <g> 1 g 1 <g» 1, 

A 2 2 = {±l/V2)a 3 g CT 3 g CT 2 <g) 1 g 1 g 1, 

A 23 = (±l/-\/2)er 3 g ct 3 g ct 3 <g) ct 1 g 1 g 1, 

A 24 = (il/v^)^ 3 0(7 3 0(7 3 ®(T 2 0l®l, 

A 31 = (il/v^CT 3 g CT 3 g CT 3 g CT 3 g cr 1 g 1, 

A 32 = (±l/\/2)cr 3 g) CT 3 ® CT 3 g) CT 3 g cr 2 ® 1, 

A 33 = (±l/\/2)cr 3 g) cr 3 g) cr 3 g) (T 3 g) CT 3 g) ct 1 , 

A 34 = (±1/\/2)ct 3 g) cr 3 g CT 3 g) CT 3 g) CT 3 g CT 2 , 


search for the ground state solution is much more com¬ 
plicated because we have to solve systems of 64 coupled 
partial differential equations for the 64 components of 
the wave function. In this extension we will also restrict 
ourselves to solve some reduced models with particular 
assumptions on the bosonic degrees of freedom. 

VI. SOLUTIONS TO THE EXTENDED MODEL 

Using the matricial representation in eq. (1^1) we need 
to solve, again, systems of matricial partial differential 
equations. In this case, every supercharge and gauge 
symmetry operator is a (64 x 64) operator matrix acting 
on a 64 component wave function IT). 

As an example we begin with the following model. We 
select two variables to be different from zero, (j>\ = x , 
t f >2 = V- In this particular model, we see from eq. ED 
and eq. (l23l) that it is not necessary to impose the re¬ 
strictions G a |T) = 0 for any G a . We need then to solve, 
for example, Qi|T) = 0. As in the previous models the 
components of the wave function can be identified and 
one gets a solution with just two independent compo¬ 
nents 0i, and 0 2 . The independent components (0i,02) 
appear finally coupled in only two differential equations 


where a l ,i = 1,2,3, are the Pauli matrices. The (4 x 4) 
Dirac representation for the T m matrices we use is 



dip 2 
1 dx 
.dip! 
1 dx 


dip2 

+ -w- Sxyipi = 0, 

dy 


diPi 

dy 


3 xyip 2 = 0. 


(25) 


These two coupled components (0i,0 2 ) are related to 
other six components of the wave function in the follow¬ 
ing manner 


In this case / is the 2x2 identity matrix. Then, for 
example, the first supercharge takes the following form 

Ql = — Al 3 7rD — A 23 7T 2 — A 33 7T3 — Al47T 2 (24) 

- A 24 7T 2 - A 34 7T 2 + An7T 3 + A 2 l7T2 

+ A 3 i7r| - iA 13 Wi - iA 23 7r 2 - zA 33 7t| 

- igAn(<f>l4>i - 0 3 0 2 ) - igA 2 i{4>\<p\ - <p\<pi) 

- igA 3 i(cp\(p2 - 4>\(p\) - gA 12 ((pl<pl - 0 3 0 2 ) 

- sA 22 (03<£i - (p\<t>l) ~ gA 32 {(p\4>l ~(p\(p\) 

+ igA 12 {cpl(p 3 - (pl$) + igA 22 ((p\(p\ - 0 2 0 3 ) 

+ igA 32 (<pl<pj - cpl(p\) + gA 13 {(j)\4>l - <p\(pl) 

+ ffA 23 (0^0f ~<P\4>1) + <?A 33 (<Ai<^ 2 - 020l) 

- igA 13 ((pl(p\ - - igA 23 (4>l4>\ - <p\(p\) 

- igA 3 3((pl(p 2 ~ <Pl<Pl) + <?Ai4(0 2 0| — 0302) 

+ “ 0i0i) + 9^3i{(pl<p 3 - 0|0?)- 

In this extension we have a richer structure and we can 
search for more interesting potentials. We notice that the 


08 = *01,07 = *02,019 = ~iip2, (26) 

020 = * 01 , '021 = - 01,022 = 02 • 

There are also other pairs of components coupled in the 
same way, here we write all these pairs of coupled com¬ 
ponents 

(01,02), (03,04), (09, 0io), (011,012), (27) 

(033,034), (035,036), (041,042), (043,0440 

If we identify all these pairs as one, then all the rest of 
the components are related to the only independent pair 
(0i,02) in a similar way as in eq. (1261) and the wave 
function has only two independent components. This, 
now exact, reduced model obtained following our gen¬ 
eral procedure would correspond to the only two compo¬ 
nent one guessed in [2C| to analyze the spectrum of the 
Hamiltonian eq. ©• Our quantization approach leads 
to similar features. It allows however an exact solution 
which satisfies also the Hamiltonian operator including 
the fermionic sector at once. A continuum spectrum of 


the Hamiltonian is an expected feature in the connec¬ 
tion between the supermembrane and DO-branes 0,m. 
For several possible choices of the bosonic variables one 
gets trivial solutions. Here we show solutions to some 
non-trivial models. 

A. Solution 1 

Another choice of a model with three bosonic compo¬ 
nents is the following 

4 >\= x, 4> 2 = y, 4>\=c = constant (28) 

with all other bosonic variables equal to zero. In this 
case, we should impose the condition G 31 H/) = 0, with the 
consequence that many components of the wave function 
vanish. The only components that are different from zero 
are 

01,02,03, 04, 021,022, 023, 024, (29) 

041, 042, '043 , 044, 061,062, 063, 064, 

and solving Qil 1 !') = 0, leads to the solution which is in 
this case 


The rest of the components are also coupled in the same 
way and give the same solution, these pairs of coupled 
components are 


(01, 

037), 

(02, 038), 

(03,039), (04, 040), 

(05, 

033), 

(06,034), 

(07,035), (08,036), 

(09, 

045), 

(010,046), 

(011,047), 

(012,048), 

(013 

, 04l), 

(014,042), 

(015,043), 

(016,044), 

(017 

,053), 

(018,054), 

(019,055), 

(020,056), 

(021 

,049), 

(022, 050 ), 

(023, 051), 

(024,052), 

(025 

, 061), 

(026,062), 

(027,063), 

(028,064), 

(029 

,057), 

(030,058), 

(031,059), 

(032,0600 


It can be seen that there are only two independent eigen¬ 
vector solutions, one of them is normalizable and the 
other one is not. In the next section we will see that 
this last solution can be interpreted as the supersymmet¬ 
ric quantum solution of a cosmological model. 

VII. CONNECTION WITH SUSY QUANTUM 
COSMOLOGY 


01 


A exp 




021 


—iAex p 



01 = 02 = 03 = 04 = 041 = 042 = 043 = 044 
021 = 022 = 023 = 024 = “061 = “062 = ^063 = 


~y 2 ) 

(30)' 

-064- 


with A a real constant. 


B. Solution 2 

The next one dimensional model has a more interesting 
form, it has four bosonic components and three of them 
are constants. The choice of variables is 

02 = x, 02 = 03 = c, 03 = d, (31) 


where c and d are constants and the rest of the bosonic 
degrees of freedom are set to zero. The conditions 
GqIT) = 0 are not required and the solution to Qil'F) = 
0 leads to pairs of coupled components, for instance 
( 01 , 037 ) that can be placed in a two dimensional base of 
two eigenvectors in the form 


01 

037 




+ C 2 


1 

i 




(32) 


As it was discussed in the introduction, M(atrix) the¬ 
ory is related with the theory underlying the five consis¬ 
tent theories of strings and 11 -dimensional supergravity, 
namely M-theory. It is then expected that the matrix 
models would be related to gravity as well. A relation 
between gravity and matrix models is that these models 
describe the properties of Schwarzschild black holes af¬ 
ter compactifications on torus of various dimensions [38l — 
l45j . It has also been proposed that one can construct 
a cosmological model corresponding to matrix theory, in 
j!4(| the authors study the relation between this matrix 
model and Newtonian cosmology. In [42$, using a ho- 
mothetic ansatz, the authors were able to get the Fried¬ 
mann equations from the classical equations of matrix 
theory. Here, it is important to remark that in all these 
cosmological models, the fermionic sector is absent and 
that they have been worked out only with classical ma¬ 
trix equations. We are dealing with the quantization of 
the matrix model of the supermembrane directly related 
to the matrix model of DO-branes. Then we can expect 
that, given that we are working with the quantization 
of a supersymmetric model, our solutions would be re¬ 
lated to the solutions that have been obtained in differ¬ 
ent approaches to SUSY quantum cosmology [27l [28l. [30l— 
[32} . As it was pointed in j28j, there are several ways 
to carry out the quantization of a supersymmetric cos¬ 
mological model. One of this options is to use the su¬ 
perfield approach, generalizing the minisuperspace of the 
FRW cosmological model (and other cosmological mod¬ 
els of interest) introducing the corresponding superpart¬ 
ner coordinates of time t —> (f, 77 , fj) and working with 
the supersymmetric generalization of the FRW action in 
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terms of the superfields corresponding to the scale factor 
and the lapse function respectively R(t) —>• M(t, 77 , fj) and 
N(t) —»• N (t, r /, fj) ®,H1. We can actually note that one 
of our solutions resembles the supersymmetric cosmolog¬ 
ical solution found in [32j with this last approach. 

From the supercharges and the SU(2) symmetry oper¬ 
ators we can calculate the Hamiltonian of every model 
we have solved using eq. ©■ For the solutions in 
the extended 4-dimensional model the Hamiltonians are 
(64 x 64) non-diagonal matrix operators whose diagonal 
terms are related with its bosonic part. Let us first take 
a look at the second solution of our extended model eq. 
(l32l) but identify now the variables and constants in this 
case as 


In the superfield formulation one introduces the su¬ 
perpartner coordinates of time in the superspace t —> 
( t,p,fj ). The superspace valued Taylor series expansion 
of the superfields is 

N (t, r), fj) = N(t) + iry4>(t) + ifjip{t) + V ( t)pfj , (38) 

R (i, 77 , fj) = R(t) + ipX(t) + ifjX(t) + B(t)pfj, 

where ip(t),ip(t) are the complex gravitino, and V(t) 
is a U( 1) gauge field. B(t) in the superfield R(t,p,fj) 
is an auxiliary degree of freedom and A(t),A(f) are the 
fermionic partners of the scale factor R(t). The super- 
symmetric generalization of eq. (13611 is 


4 = m, = = = (33) 

where k is a constant, c is the speed of light and G is 
the gravitational constant. We will also recover for this 
solution the constant h. The bosonic (diagonal) part of 
the Hamiltonian is 


ir 2 kc° 2 2Mc 4 i/k„ , .0 o . , „ 

Ha = p R + ~q2^ q - R M c ) * = 1, 2,..., 64 


(34) 

and the solution in eq. (l32l) for the coupled (ipi,ip 37 ) 
components of the wave function become 


^37 


= Ci 


c 2 


1 

— i 

1 

i 


exp 


exp 


_ ( Mc r? V^ c3 

V h R 2 gh r 


Me . 


-R- 


^ D 2 


2 Gh 


- ir 


(35) 


Now we review how this solution arised in the con¬ 
text of SUSY quantum cosmology. Following [32, [Hj we 
rewrite here the supersymmetric action of the FRW cos¬ 
mological model with a perfect fluid of barotropic equa¬ 
tion of state p = 7 / 9 . For the case of dust 7 = 0 and null 
cosmological constant A = 0 such action is given by 


S = 


2NG R {™) +» K 4 R ~ NMc2 


dt (36) 


where k = —1,1,0 is the curvature constant. N(t) and 
R(t) are the lapse function and the scale factor respec¬ 
tively and M is the mass parameter of dust. This action 
as given is invariant under reparametrizations of time 
[t —> t + a(t)] when N(t),R(t) transform in the following 
way 


r D dR 
on =a——, 
dt 


SN = 


d(aN) 

dt 


(37) 


S = 



RZl^RZl^R 



- McR 


dpdpdt, 


(39) 

where D n = --§= - irp^ and D v = ^ are the 

supercovariant derivatives of the global “small” super- 
symmetry of the generalized parameter corresponding to 
the time variable t. When the canonical formalism is 
applied to the supersymmetric Lagrangian it was shown 
[3lL 1321 that the classical total Hamiltonian is given by 


Hc = NH+ l -^S- l -yS + \vF, (40) 

where H is the Hamiltonian of the system, S and S are 
the supercharges and F is a U(l) rotation generator. N , 
if), ip and V are now Lagrange multipliers. H, S, S and F 
are first class constraints, hence the physical state wave 
functions *F (R) are obtained from the conditions 

H'&(R)= 0, SV(R), S^(R) = 0, F*(R)= 0. (41) 

These constraints have a dependence on R(t) and 
A(t),A(f). The explicit dependence of the Hamiltonian 
and supercharges on this variables is the following 


„ G 2 nc 4 R 
2 c 2 R nR 2 G 


M 2 G 
2 R 


+ c 2 M^/k 



MG 
2 cR 2 


XX, 


s = 


iG 1 ' 2 c 2 ^R x ' 2 

cR 1 / 2 7rR G 1 / 2 


MG 1 ' 2 \ 

A 


s = 


iG 1 / 2 c 2 v ^i? 1 /2 

cR 1 / 2 ' nR+ G 1 / 2 


MG 1 / 2 \ - 

) A ’ 


(42) 


where ttr is the canonical momentum of the variable 
R. As a consequence of the superalgebra between the 
Hamiltonian and the supercharges, the physical states 
are found by applying only the supercharge operators to 
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the state wave function T(i?). Using the following matrix 
representation for the fermionic variables 

A = Vhcr_, A = -Vha + , a± = ^(cr 1 ± a 2 ), (43) 

T(i?) becomes a two component wave function T(I?) = 
(T 1; T 2 ). The corresponding superquantum solutions 
1321 for these two components are 


Ti 

T 2 


Cexp 
(7 exp 






(44) 


We can see that our solution eq. (1551) and this one ob¬ 
tained in the superfield approach eq. (1441) . are indeed 
the same. The diagonal part of the matrix Hamiltonian 
for every matrix model is associated to the bosonic sec¬ 
tor of the system. If we want to find the solution to 
Hty(R) = 0, we can rewrite the Hamiltonian constraint 
in eq. (l42l) and see that the bosonic part of this Hamil¬ 
tonian, which is the one that does not depend on the 
fermionic variables (A, A) is the same as the diagonal part 
of the matrix Hamiltonian in eq. (1341) of our particular 
matrix model. We have then found the same solution in 
SUSY quantum cosmology and in matrix theory follow¬ 
ing our general procedure to solve matrix theory. 

The behavior of the Ti component is the expected 
when R —> R sup , where R SU p is the maximum radius 
of the universe. One of the components converges as 
R —> oo while the other one diverges, but the convergent 
solution can be isolated with the right choice of the con¬ 
stants. The scalar product of the Ti component in the 
measure R 1 / 2 dR is normalizable 


1 = C 2 f ^ ^lR 1/2 dR, (45) 

Jo 

= C 2 J q RSUP (e 2MR ~^) R^dR, 

and in both cases the norm of the wave function can be 
defined with the inner product 

(TIT) = [ P T t T R 1/2 dR. (46) 

Jo 

This review allowed us to know how to get a supersym¬ 
metric quantum solution for a particular cosmological 
model in the superfield approach, and we were able to 
relate it with ground state solutions we have got in the 
matrix model, in particular for the second solution of the 
extended model. It was then important to show in some 
detail how these supersymmetric quantum solutions are 
obtained and how they are related with possibly normal¬ 
izable solutions in matrix theory. We want to remark, 


after we showed the coincidence of the solutions, that 
this makes evident certain connection between solutions 
in matrix theory and solutions in cosmological models 
arising from supergravity. 


VIII. CONCLUSIONS 

We have considered finite N matrix models, the rel¬ 
evant operators are the supercharges and those corre¬ 
sponding to the SU(IV) symmetry. The Hamiltonian op¬ 
erator is a consequence of them. The fermionic degrees 
of freedom are represented as Dirac-like gamma matrices 
[27l l28j . By means of this, we were able to propose a 
general procedure to find exact ground state solutions to 
SU(2V) matrix models. We have shown here some super- 
symmetric reduced models and found their exact ground 
state solutions being some of them normalizable. The 
method we propose allows the possibility to find exact so¬ 
lutions by means of the constraints related to the SU(-ZV) 
invariance and the corresponding supercharges, contain¬ 
ing both the bosonic and the fermionic sectors. Even 
though we have solved reduced models, the extension to 
any dimension and to any SU(IV) gauge invariant group 
is possible following the proposed procedure. One could 
then search for exact solutions to the full 9-dimensional 
matrix model and for any finite N although powerful 
computational tools would be needed because the dimen¬ 
sionality of the matrices involved make the problem hard 
to handle. Given the structure of the supercharges for all 
the allowed dimensions as functions of the momenta and 
the potentials of the variables involved, we expect that 
some solutions to the lower dimensional models would 
also be solutions to certain higher dimensional ones. This 
is the case for the solutions we have shown here. Some 
of the solutions found to the 2-dimensional model are 
also solutions to the 4-dimensional extension. A reduced 
model guessed and investigated in 0 results as an ex¬ 
act sub-model in our formalism. A particular interesting 
result is a relation between one of our solutions in the 
extended model with supersymmetric quantum cosmol¬ 
ogy. This guarantees that the reduced matrix model is 
supersymmetric. These relations reinforce the method of 
quantization we propose in this work for finite N matrix 
theory 0- Besides, we encounter possible relations be¬ 
tween several topics connected with our work; on the one 
hand, according to our results, there is a connection be¬ 
tween reduced matrix theory models and SUSY quantum 
cosmology, also, already in [47] the FRW cosmological 
models have been related with classical bosonic matrix 
theory. Even Schwarzschild black holes [38j - l4l| and their 
thermodynamical properties [T3 [d~5l | have been obtained 
by compactifications in matrix theory. On the other 
hand, in a recent work 0] the authors show evidence 
of the conjecture relating supergravity and the dynamics 
of a spinning particle moving in an infinite coset space by 
means of the quantum dynamics of the supersymmetric 
Bianchi IX cosmological model [0 and the operational 
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structure of the constraints revealed a hidden hyperbolic 
Kac-Moody structure. It would be of interest to be able 
to identify reduced matrix model solutions associated to 
the Bianchi SUSY quantum cosmological models, in par¬ 
ticular the Bianchi IX [28|,|4H , and by means of this try to 
be able to find a relation between the Kac-Moody struc¬ 
ture and matrix theory. An analysis of our work regard¬ 
ing several of these possible connections is the matter for 
future research and is beyond the scope of this work. 
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